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FOREWORD 
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The research reported in this document is concerned with the design of 
geophysical fluid dynamics (GFD) experiments of major scientific 
importance which require the unique "low-gravity" environment of Earth 
orbit to be successfully performed. The experiments offer the opportunity 
to study baroclinic fluid flows on rotating spherical surfaces with an 
imposed radially directed body force to simulate a gravitational force 
field. Comparable terrestrial baroclinic GFD experiments have been 
restricted to cylindrical geometry, and thus to the study of g-plane 
flow situations. Accordingly, Ekrth orbit offers the potential for a 
major step forward in GFD research which in turn may lead to new under- 
standing of the planetary-scale dynamics of atmospheric flows. 
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1. Introduction 


Meteorologists and astrophysicists interested in large-scale pla- 
netary and solar circulations have come to recognize the importance of 
rotation and stratification in determining the character of these flows. 
In particular, the effect of latitude-dependent Coriolis force on non- 
linear convection is thought to play a crucial role in such phenomena 
as differential rotation on the sun, cloud band orientation on Jupiter, 
and the generation of magnetic fields in thermally driven dynamos. Most 
theoretical and all experimental work on these problems has in the past 
treated only local curvature effects — the mid- latitude or equatorial 
3-planes of meteorology being well-known examples. In fact, terrestrial 
laboratory experiments have only been able to study 3-plane flows in 
situations where the stratification is extremely simple, usually two 
layers of fluid of slightly different density (see Hart, 1972). Even 
here only local curvature effects occurring near the poles can be 
modelled. One cannot model equatorial regions where the rotation axis 
is perpendicular to gravity because then the earth's gravity field 
becomes an oscillatory force in the coordinates rotating with the body. 
Other body forces proportional to density can be generated in the lab- 
oratory. One simple example is the centrifugal buoyancy force which 
is perpendicular to the rotation axis. Thus if a body is rotated rapidly 
enough, this centrifugal force can dominate regular gravitational effects 
and a dynamical approximation to the equatorial plane can be achieved. 
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This effect has been used by Busse, et. a1 ., (1974) to study rotational 
effects on equatorial convection. However, this analogy extends to non- 
equatorial latitudes only when the body rotates so fast that rotational 
constraints are dominant everywhere (e.g., as the Taylor number T. , defined 
below, becomes “). The values of T, appropriate to the astrophysical 

a 

2 6 

flows we wish to study are thought to be in the range 10 to 10 (Gilman, 
1974). The latter is a strong (possibly but not dominant) rotational 
limit, while the former is a weak rotation case. Clearly we need to study 
both limits in order to gain understanding of geophysically relevant flows. 
The key to physical experimentation on this subject is the generation of 
buoyancy forces which are radially rather than axially or uniformly di- 
rected. This report describes the theoretical design of such an experi- 
ment along with results of preliminary laboratory tests of a prototype. 

2 . Theoretical foundations 

In trying to construct a physical laboratory analog to a geophysi- 
cal flow it almost always proves impossible to maintain geometrical and 
dynamical similarity in all respects with known materials and techniques. 
The case of deep, convecting, rotating, planetary or solar atmospheres 
is no exception. We take the point of view that we will make a labora- 
tory simulation of a system which is in many respects similar to the 
natural ones we are interested in, but in which already some processes 
are neglected. The fundamental system we take is that of Boussinesq 
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liquid contained in a rotating spherical annulus with thermal forcing 
(which may be latitude dependent) imposed at inner and outer spherical 
surfaces. The configuration is shown in Figure I. The body rotates 
at rate S2 and contains a heated annulus of k, v, a-type fluid. The 
following parameters describe the system and fluid. 


r 

0 

<t> 

Q 

K 

V 

OL 

«i 

«o 

p 

H. 

Po 

T 


Table I 

Fluid parameters and geometry 
radial coordinate 
longitude 
latitude 

basic rotation rate 

thermal diffusivity 

kinematic viscosity 

coefficient of volume expansion 

inner radius 

outer radius 

pressure 

velocity vector 

background mean density 

temperature relative to ambient 
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Figure 1. Geometry of the rotating spherical shell experiment. The 
outer shell at radius is maintained at temperature T^, voltage T = 0. 
The inner shell at radius is at temperature T^. voltage V sin cot. The 
applied frequency is 60 hz. The temperature difference AT = T^. - T^ 

(> 0 for unstable convection). The whole apparatus is rotated about 
its axis. The interior is fed via a supporting neck. 
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The dimensional equations In this Boussinesq fluid model are: 

^ <^Cr) o( ^ V 

dJ = k T 

cit ^ 

= O. (2.3) 

Note that the buoyancy force acts In the radial direction with a given 
gravitational acceleration g(r). We wish to construct a laboratory 
analog of the flow situation described by these equations. 

Let us first note the well known assumptions used In deriving 
(2.1) - (2.3) from the full non-linear Navler-Stokes equations. 

1. The Mach number Is small 

2. The height scale for the motion Is small compared with the 
atmospheric scale height. 

In addition to applying these Incompressible equations to a deep compres- 
sible gas we also assume that small scale turbulence can be represented 
by uniform "eddy" viscosities and diffusi vltles. One might want to 
add small scale turbulence (occurring at scales « d), magnetic fields, 
radiative transfer, and compressibility, to make this mathematical model 
more realistic. However the possibility of constructing a laboratory 

•ff 

For conciseness we denote frictional terms as "V^" recalling that In 
different coordinates the divergence of the viscous stress tensor In- 
cludes new terms. 


( 2 . 1 ) 

( 2 . 2 ) 
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analog with all these added physical elements is extremely remote. 

Equations (2.1) - (2.3) apply to subsonic convective motions of 
a shallow rotating atmosphere when a is replaced by the inverse 
background potential temperature and T by 0. The equations stand as 
written for all subsonic oceanic motions except where different dif- 
fusivities of salt and heat are important. Thus a laboratory model 
based on these equations will have a wide range of applicability to 
many different geophysical situations. 

Another reason for starting with eqns. (2.1) - (2.3) is that 
they serve as the governing equations for most theoretical studies of 
planetary and solar circulations and contain much of the key physics 
without added complexity of the other effects. While it is admittedly 
risky to extrapolate results derived from these Boussinesq equations 
to atmospheres, they do contain descriptions of many of the processes 
believed to be of fundamental importance. Scale interactions via non- 
linearity, convective driving, and the latitude dependent Coriolis 
forces are all present. That we now have only a meagre understanding 
of thermally driven flows described by (2.1) - (2.3) suggests that a 
laboratory analog of these somewhat simplified equations should be 
very useful in verifying existing theoretical results and in providing 
useful scale information for generating further approximations in the 
solutions of these equations. 



It is not altogether obvious that a laboratory analog of (2.1) - 
(2.3) can be made. The principal difficulty arises in generating the 
radial buoyancy force. We propose to do this in a system which takes 
advantage of the temperature dependence of the dielectric constant of 
certain liquids. Consider a dielectric liquid contained in the spherical 
annulus of Fig. 1. The dielectric constant is 

e = [l - !i' (T-T„ )] , , 


where is the ambient dielectric permittivity, T^ the ambient temp- 
erature (about 22“C in the lab) and y the thermal coefficient of the 
permittivity. The equations governing this dielectric fluid flow in an 
imposed electric field E(^^) is, in the laboratory, 

e^'f ^ * j ^ / = 

•2 

- M. V e 

^ ’ (2.5) 


^ + (u . 7)T - ^ ^ 


e«cp 


HT{y'7)E 

a — >> 

eo Cp 


( 2 . 6 ) 


(2.7) 


V ‘ U - O 


7 
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where Cp is the heat capacity of the dielectric liquid and z the direc- 
tion of the earth's gravitational acceleration. $ represents viscous 
heating which is neglected. Now note that if the dielectric heating 
terms are small enough, and if we can construct a situation where 

BV€ ~T V -f(h) ^ 

then the laboratory model will be an analog of the Boussinesq equations 
on a rotating sphere with radial gravity. 

Table II lists properties of a typical silicone oil. 


V 

K 

a 


P 
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0 


c 
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Table II 

Properties of Dow Corning 200 fluid 


.ol cm /sec 


6.4 X 10"^ cm^/sec 
1.17 X 10"^ °C"^ 

3 

.873 grm/cm 
;.0 X 10“^ °C"^ 

.45 cal/grm °C 


10"® m^/sec 
6.4 X 10"® m^/sec 
1.17 X 10"® °C"^ 
873 kg/m® 

1.0 X 10"® °C"^ 
1.9 X 10® j/kg °C 


This fluid has a Prandtl number P^ = v/k = 15. 

For various reasons we will be setting up an ac electric 
field between the conducting boundaries of the spherical annulus. 
The frequency CO of the field will be constant and much greater 
than that of the thermally driven fluid motions whose time scale 
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will be on the order of the rotation period or the advective time 
scale, whichever is smaller. The motion frequency will be less than 
1/sec. Over these time scales of 1 second or larger the first thermoelectric 
term in (2.6) has zero average and can be neglected. The second has 
magnitude “T d f iJ 1/ ^ 

V5 eCpCT;--rO ^^ d ' V 

relative to the conduction term. For an applied temperature difference 

fi — i 

of 2 C, a field of 10 V/m, and a 1 cm gap, this ratio is 10 or less. 

If the dielectric fluid is lossy an internal heat generation term 

of form 00 E CO Li should be added to the left hand 

side of (2.6). For silicone oil the loss factor L is so low that this 

- 5 

term is of order tO with respect to the conduction term. The fluid 
resistivity is so high as well (10^® ohm/ cm) that the Joule heating is 
infinitesimal. Thus to an excellent approximation the thermodynamics in 
this laboratory situation ore essentially identical to that in the original 
Bousinesq model where only diffusion, advection, and time changes of 
temperature were considered. 

Now let us look at the momentum equations. We have neglected 

free charge forces cr^E. li^ the driving field frequence w is much 

-1 ”2 

greater than the relaxation frequency ~ (e ■>- resistivity) - 10 sec 

this force is negligible. This is why we must drive with an ac field 

with frequency ~ 60 hz. Consider the forcing terms we are retaining, 

namely . 

^ Vp -h <^eo<^T 72 ^ E ^ VT /a, 

appearing on the right hand side of (2.5). In a divergence-free fluid 
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with no-slip boundary conditions the dynamics are determined by vor- 
ticity equations obtained by taking the curl of (2.5) and (2.1). The 
vorticity equations will be identical except for differences in curl 
We have 

CijrJ £ = 

' ( 2 . 8 ) 

2 2 

Thus if E = E (r) and the background gravitational acceleration is 

small, the system of equations (2.5) is equivalent to (2.1) where the 

original p.g(r)aVr in (2.1) is replaced by VE^(r). 

° 2 
We can look at this another way. Since the boundary conditions 

don't involve the pressure we can rewrite a new composite pressure 

p p _ E^Or)^ T /X 

Now 

■f = - + g^„o(T7z - VE 0-) 

3L 

= - 7p„ T(<i^ + 

where is the laboratory geopotential 

and (()g the equivalent electromagnetic geopotential, 

- ^ /x . 

Thus the primitive equations are equivalent since on an astronomical 
body (J) = apQG(r). 
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Experi mentally we impose an alternating electric potential 
between inner and outer shells. We can ground the inner shell and 
apply if) = V sinwt to the outer shell. The average potential if; in the 
liquid is then 



Using this we can compute the ratio of the terrestrial buoyancy 
force f^ to the electrodynamical buoyancy force f^^. We have 


(2.9) 


This equation shows the need to make Rq-R^ and R^ itself small. 

We would like to make very large so that the buoyancy force will 

be effectively radial. For = 2 cm, Rq = 3 cm, 

= .002 (^kv) / ci 

using fluid data from Table II. Clearly large fields 0(50kv) are neces 

sary if is to be large in the terrestrial laboratory where 

2 

g ~ 10 m/sec . Unfortunately dielectric fluids tend to be electrically 
unstable at these fields. However in an orbiting laboratory with 
g -.001 m/sec a 2.5 kv field gives a f^/f^ ratio of 12. 

Further discussion isfacilitated if we nondimensionalize our equa 
tions using the scales </d for velocity (d = Rq"R-{)» d /< for time, 

AT for temperature (AT might be the difference between inner and outer 
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2 2 

applied wall temperature) and k /d for pressure. These scales are sub- 
stituted into (2.5) - (2.7) after the thermoelectric terms in (2.6) and 
the centrifugal terms in (2.5) have been dropped. This latter action 
is consistent for laboratory flows where 

1 

^ 2 

With g ~ 1 m/sec (the electrodynamic "accelerator"), - .01 m, this 

condition implies that n « 10 rad/sec. We will be concerned with models 

where n ~ 1 so will be small, 
r 

We define the characteristic dimensionless parameters of this 

problem: o 

g oATd'^ 

R_ = Rayleigh number = 

o KV 


> o 


= Prandtl number = v/k 


( 2 . 10 ) 

( 2 . 11 ) 


= Taylor number 

a 


4n^d^ 


V 


3 =Ri/d 


( 2 . 12 ) 

(2.13) 


If IT is the electrodynamically modified dimensionless pressure, the 
governing equations are: 

j. + PT-''^{^xu ) 

'ha \ 

= _ V TT ( 7 V -i- ) 


4 - 


A>a PrTi/7C ^ P.R-T P 


(2.14) 
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= v^“r (2,15) 

ci± 

v*^ = o , (2.16) 

In the above equations is the electrodynamic acceleration 

. = +;S)^~,oo3 

and 

X = 3°/^ ' ^ 

In the terrestrial laboratory x —02 at 8 kv, in orbiting laboratory 
X - 20 at the same voltage. Thus we see the great advantage (or even 
necessity) of doing these experiments in space. Suppose the boundary 
conditions at r = B» 3 + 1 (inside and outside shells) are 

u = V = w = 0 
T(3) = 3+1 

T(B + 1) = 3. (note AT = T(3 + 1 ) - T(3) = -1) 

Then T = 43(3 + l)/r, u = o, it « solve the basic equations when 

X -»• <». This basic state of no motion and radial (unstable) temperature 
gradient can be superimposed onto the field variables. In spherical 
polar coordinates the equations with u, v, w being the east (0), north 
(+4>) and radial velocities respectively are; 

^ + U . 7U -4- UuT ^ u i/ faff ^ ^21 

dt r h ' f-coS(^ dO 

^ /2. 

- ^ COS ^ V S)/?^ j 

It 

Numerical values computed with R^=2cm, d=lcm, and fluid data from Table II. 
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^ ^ s //I ^ 5 TJ' ^ 

P~coT^ BiP A*C<3S*<Z< To ^ 7\ap^ 


^ ^ a ‘ vv-^ T^ u ^-h^d - _ _L aZT P IZ \ • J,, 

■ n ^ ^ :?co5^ ji:/ ^ 

P^T COS^/tC + /^Py.^(^+l)co5^/K-X- ^ 

iktr 6/ . _ Z/*" _ _ 97T P J- 

-b6 7^ ~ - iT * cc?s^6^ 

-+ P 5 17 6^ - - -^- Bl/cosd _ a ^ T 

^ r cc>^ ^ /> *• sih (^ '^0 J (2.19) 

+ (Hi/ /?.PrT -I- R,Pr s'/n (j> ~T /pC f Rn?r ^ 


( 2 . 18 ) 




U • V T = -H V^T ^ 


( 2 . 20 ) 


j_ \r^uA t 


cos ^ 9^ 


^ T^COs 4 -h 


_! 2i! = o 

fco&d * 


C2.21) 


where 



zr ^ . Ud 

/“ tcos^ 


U . 7" 
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3. Conwients on the electrohydrodynamics 

In the previous section we showed that for large values of x the 
thermally driven motion of a dielectric fluid in a spherical capacitor is 
described by the same equations as that of a thermally driven rotating 
shallow atmosphere or reasonably deep ocean on a spherical planet, when 
the small scale motions for the latter are modelled by constant eddy 
viscosity. Several assumptions were made in the representation of the 
electrohydrodynamical effects in the fluid which were important in ob- 
taining the desired result that a radial gravity field could be simulated. 
Although the subject of electrohydrodynamics is a relatively new field 
there have been several previous studies which lend support to the assump- 
tions used here. 

As mentioned above, there are several surface instabilities related 
to surface electrical stresses on a dielectric liquid with a free surface. 
These are discussed in the review by Melcher and Taylor (1969). Suffice 
it to say that they have no geophysical analogs and are undesirable in 
our model experiment. They are eliminated by applying rigid electrically 
conducting boundary conditions to the fluid. Since this is the natural 
configuration for our spherical capacitor we will have no trouble with 
the surface modes of instability. 

Among the bulk fluid instabilities possible when the dielectric fluid 
is subject to both thermal and electric fields are (a) the electrohydro- 
dynamic Rayleigh-Taylor instability that we desire in our study of solar 
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and planetary convection and (b) an instability related to free-charge 
effects Turnbull (1971a). We have assumed that instability (b) can be 
eliminated by applying an ac field of high enough frequency. Turnbull 
and Melcher (1969) and Turnbull (1970) have shown by comparing theoreti- 
cal and experimental results that the Rayleigh-Taylor instability 
in a plane layer of fluid between horizontal conducting surfaces be- 
haves as described by the Boussinesq equations neglecting thermo- 
electric and free-charge effects. In their experiments the gradients 
of E are small and related to the weak temperature dependence of e 
(e.g., V • 5 = 0 yields an E with gradients of order yAT). In our 
experiments a strong E gradient is imposed geometrically so that free 
charge effects will be relatively weaker in comparison with polariza- 
tion force effects. Their experiments were done using a silicone fluid 
quite similar to the one we propose to use and the good agreement they 
achieve encourages us that we are properly representing the electrohydro- 
dynamics. A similar study of the effect of a non-uniform alternating 
electric field on a thermal boundary layer was reported by Turnbull (1971b) 
and this too gives us confidence in the reduced equations (2.14) - (2.16). 
Further supportive evidence comes from the study of Chandra and Smylie 
(1972) who show that the convective instability of a differentially 
heated dielectric fluid in an alternating axially symmetric field occurs 
at the predicted critical Rayleigh number computed using equations like 
(2.14) - (2.16) with x “ as a starting point. 
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Based on these studies we have reasonable confidence that the elec- 
trohydrodynamical behavior of the silicone fluid at moderate temperature 
differences and large 60 hz fields is described by our model equations. 

It remains to verify this for our particular geometry. This is largely 
a question of how to get adequate data. We discuss this in following 
sections. Note however that we have also assumed that £ is a function 
of r only. In a perfect spherical capacitor this is obvious, but in a 
real experiment we need to have a feedthrough tube to maintain the inner 
sphere at constant voltage and temperature. Near the "neck" then 
^ ^ (r, <|)) and the term VE x 7T will not be zero for the basic state 

T (r) = -a(l+B)r as it should be. Thus we must conduct our experiment 
far enough away from the neck so that E is primarily radial. This is 
accomplished by performing experiments only in the upper hemisphere 
opposite the feedthrough. The lower hemisphere is isolated by a radial 
plastic disc at the equator. For small but reasonable neck sizes and 
B ~ 3 the field, as calculated by numerical solution of Laplace's equa- 
tion, will be about 95% radial in the entire upper hemisphere. Note 
that the lower voltages required when the experiment is run in a low g 
environment means that the neck can be substantially smaller than that 
needed for insulation in the terrestrial laboratory. Thus it is probable 
that in the zero g environment more than a hemisphere could be studied. 
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4. Geophysically relevant parameter ranges 

To discuss what values of R_, T^, 3 and would be most interesting 

a a I 

to use in the experiments requires an appreciation of what astro-geo- 
physical problems are important. A complete review of what is known 
about convection in a rotating spherical shell and of what we would like 
to know is beyond the scope of this report. Let us briefly look at the 
conceptually simplest case of uniformly heated smooth surfaces to illus- 
trate the kinds of questions we can hope to answer with this experiment. 
Obviously, more complicated thermal boundary conditions and surface irregu- 
larities could be included and studied experimentally. The types of prob- 
lems are perhaps obvious to meteorologists and oceanographers^ The ques- 
tions arising in the study of non-linear mechanics are similar in a wide 
manifold of problems so we will focus on the simplest forced problem- 
one whichis by no means well understood. 

The uniformly, unstably heated, rotating spherical shell of fluid 
serves as a model problem for large scale planetary and solar convection. 
The principle motivation for this model comes from observations of large 
scale giant convective cells on the sun and (possibly convective) axi- 
symmetric rings at low latitudes and cells at high latitude on Jupiter. 

The study of this problem has by and large followed that of convection 
in a plane non- rota ting fluid layer. In particular we would like to 
understand the mechanism for planform selection. In the plane layer the 
preference for two-dimensional rolls or three-dimensional hexagons is 
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related to calculable non-linear stability properties. In the spherical 
rotating shell the calculations are by no means as simple partly because 
with the latitude variation of Coriolis force the global linear problem 
is non-separable in the horizontal coordinate. Under different external 
constraints (rotation, viscosity, etc.) it is theoretically possible for 
the linear convective instability to consist of banana cells (Busse, 1970) 
extending from pole to pole, axi symmetric Z)nes or belts (Williams and 
Robinson, 1973) or various forms in between (Gilman, 1972). In all these 
theoretical studies, convection first occurs as weak motions, superim- 
posed in the background state of unstable temperature gradient and no 
motion when exceeds some critical value R,^ (T^,S, P„). In plane 
non-rotating convection with rigid conducting boundaries R^^ » 1708. 

aC 

Gilman (1972) has determined R^^ for the equatorial annulus. The compu- 

qC 

tatio.ns, which are quite time consuming, show that R^ increases monotoni- 

a 

cally with T, from about 700 at = 0 to 10^ at T, = 10^. Clearly, 

a a a 

since we want to study convection, our experimental values of R^ must 

a 

exceed these numbers. 

The non-linear processes are of most interest since R. is very much 

a 

greater than R^_ for solar and planetary atmospheres. The appropriate 

aC 

values of R^ for the astronomical cases are not easy to estimate because 

a 

of poor knowledge of the eddy viscosity v and the depth d. Clearly the 
R^ appropriate for the large scale convection is very large. For the 

a 

sun R^ ~ lO^R^^, using d = lO^^cm, v = 3 x 10^^, k ~ v, and g^^/T_ = 10^. 

a dC 0 
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This represents little more than a best guess. It is clear however 
that the cases of most interest in astrophysical problems are those 
for which is many times critical. These values are easily attained 

a 

in the experiment. 

As R, increases above R^_ several interesting questions emerge. 

How long do the linear modes persist? Do the orientations of planforms 
suggested by linear theory persist, or do they become unstable and get 
replaced by a new planform? This is of particular importance in Williams 
and Robinson's (1973) model of Jupiter's cloud bands since they got axi- 
symmetric bands by setting = o in their equations I They did make 
a numerical stability analysis of their axi symmetric solutions but due 
to the enormous computation time required this aspect of the problem is 
by no means solved. In our experiments we should be able to observe 
the preferred cells' shapes and orientations as a function of R, - R^^, 

Tg, and 3 . 

As the convective flow field becomes moderately non-linear (R. - R^^)/ 

ol etc 

R _ ^ 1, the "simple" modes of linear theory will interact with themselves 
and other modes to produce a spectrum of disturbances and perhaps to 
generate a mean axi symmetric flow. The question of what kind of mean 
flow the convection cells can produce is of great importance in view of 
the observed equatorial accelerations on the sun and the giant planets. 
Numerical models (Gilman, 1973) for relatively small (R^ - R^^)/R,^ 
have shown the tendency to generate mean flows. Being three-dimensional 
and time dependent such models are complicated and cannot handle too 
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wide a spectrum of disturbance (hence they are confined to moderate 
non-linearity). The laboratory experiment should be able to verify 
these moderately non-linear results and provide information on wave- 
wave interaction which will hopefully provide data for realistic theo- 
retical approximation. In addition measurements should be made of 
mean-flows generated at strongly non-linear conditions (large 
which will probably remain inaccessible to numerical experimentation 
for some time. 

In order to ahieve the above goals we must be able to obtain key 
measurable quantities* a problem discussed below, and have access to 
the proper parameter range. For a typical experimental set up in a 

g 

zero g environment we find that R^^ can be set from 0 to 10 and from 

5 

0 to about 4 X 10 . Thus one can cover the weakly and strongly unstable 
limits for both the weak and strong rotation cases. Figure 2 shows a 
parameter plot of the regions accessible by this experiment on a diagram 
which also shows the numerically calculated neutral stability curve 

for large scale convective instability of rotating stratified flow on a 
sphere (from Gilman, 1972). The experiment can give data on flows which 
are more strongly non-linear than can currently be handled on a computer 
or by asymptotic methods. Note also that viscous effects will be con- 
fined to lengths of order d/T^^^^ or d/R^^^^ so that in the large 

a a 

(R_, T.) regions viscosity will be apparent only in narrow boundary 

ol a 

1 ayers . 



Figure 2. Typical experimental range for direct electroconvection In a spherical 
shell with uniformly heated surfaces. The solid curves show the theoretical linear 
stability curves. The hatched line encloses the parameter values accessible experimentally. 



10 



This is only one example. In general we expect to be able to study 
weak and strong rotational effects in both weakly and strongly non-linear 
flows. The Prandtl number of the silicone fluids is around 10. Ideally 
we would I’ike to be able to vary this towards smaller values. This may 
be possible with other types of fluids. 

We view this electro-convection apparatus as a tool for the experimental 
study of a large number of geophysical flows in which sphericity, stratifi- 
cation, and rotation are all important. There are a wide range of possi- 
bilities depending on the kinds of forcing imposed and on the detailed 
geometry. For example, mechanical driving and the presence of meridional 
barriers would naturally lead to flows of interest in dynamical oceano- 
graphy when the stratification is basically stable. Although we have 
focused on one particular problem, rotating convection, it is obvious 
that an experimental module in which gravity is effectively radial should 
have wide applicability. 
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5. A simple demonstration of electroconvection 

We performed a very simple experiment to demonstrate that polar- 
ization forces in a stably stratified fluid can overcome the terrestrial 
gravity. The geometry for this experiment is shown in figure 3. Two 
horizontal plates at z = 0 and z = a are maintained at a temperature 
difference AT. Thus a stable stratification is set up in which there 
is no motion . Thus any fluid motion or distortion of the stratification 
away from the form 



must come from the action of electrodynamic forces alone. This is in 
contrast to the spherical shell experiments discussed in chapter 6 where 
there is vorticity generation by buoyancy forces related to earth gravity. 
Here 

VT X V(|) ^ = 0 

A thin horizontal rod at z = z^ is maintained at a 60 Hz RMS 
voltage V= AV. The rod is thin enough so that it causes only a minor 
distortion of the background temperature field. Visualization is achieved 
by looking through the fluid layer (typically a = 1 cm, and the layer 
is 15 cm wide and deep) at a Moird fringe pattern oriented vertically. 
Thermal gradients in the horizontal direction will cause refraction 
to occur and will cause the apparent position of the fringes to move. 

For large AT = 10“C this is a rather sensitive technique and has the 
advantage that no particles need be placed in the fluid. Particles tend 
to pick up space charge at large AV and dance around. They also tend 
to settle out under g forces. Thus visualizing the thermal patterns in 
this way is preferable to particle tracking. 
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Figure 3. Geometry for the rod experiment. A statically stable thermal 
stratification is set up by heating the top lid. This is then destabilized 
by applying a large voltage to a horizontal rod located parallel to the 
y axis. 
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Figures 4-7 show the evolution of the fluid state as AV is increased. 
For AV = 0 there is only a weak thermal perturbation near the rod. This 
increases as AV is increased. At first the perturbation is primarily 
thermal. There is a very weak circulation within the elliptical thermal 
perturbation region of figures 5-6, because weak buoyancy layers are 
needed near the rod to satisfy the temperature boundary condition there. 

As AV is increased the flow breaks down into unstable plumes near the rod 
(see Fig. 7). The vertical motions are vigorous and the plumes radiate 
gravity waves into the stably stratified fluid at large distances from 
the rod. The upward motion consists of cold fluid rising, and hence is 
opposite to that generated by the g“T term in the vertical equation of 
motion. Electrohydrodynamic buoyancy forces have overcome the restoring 
force of terrestrial gravity, 'g' has been reversed. 

A formal analysis of this problem is non-trivial. We here indicate 
a rough theory for the critical value of AV needed to generate convection. 
The potential due to the rod is 




AV In. i ^1" ~ i 


rrnS 


1 


n 




(4.1) 


The diameter of the rod is p « a. It can be seen that the refractive 
patterns of Fig. 5 and 6 are strongest near the region of space where 
<j>rms is biggest, as we would expect. 

Now near the rod 
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Figure 4. Photograph taken looking through the fluid layer at a series 
of vertical lines painted on a translucent plastic end wall (at y = + A say, 
A»a). The lines 'bend' in response to index of refraction changes in 
the fluid. Here AV = 0 and there is a slight bending near the rod because 
it is a finite conductor and is not consistent with the static temperature 
distribution T = ^ z. Here a = 1 cm, z^ = 1/2 cm and AT = 20° C. The 
fluid is a low viscosity silicone oil. 



Figure 5. Same as 4, except AV = 2 kv. Note how polarization forces have acted to 
distort the thermal structure near the rod. There is a very weak recirculating motion 
near the rod. 



Figure 6. Same as 4, except AV = 4 kv. Note increased distortion as polarization 
forces are strong. Flow is still steady. 











Figure 7. Seme as Figure 4, except AV = 8 kv. / 

throughout the fluid, with cold fluid rising and hot fluid sinking. 
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where r is the radius measured from the rod center, p < r « a. Because 
of the static stability the motions will be weak (prior to electrocon- 
vective instability). If 'P is the streamf unction for motion in the x-z 
plane (w = u = + '1' ), the equations of steady motion can be com- 

bined into the vorticity equation 

o-J(T^E^) = 0 (4.2) 


where T is the thermal departure (dimensionless) from a. background .tempera- 
Ta*1ent> 


AV^ y 6, A I 


cr - . 2. ( Sirx. TT 2^0 ) 


(n 


ir 


CL 


\ ■ 2- 


and 




<X at cl 




R is the normal Rayleigh number based on terrestrial gravity. 

dl 

Now in this experiment both R^ and o are very large. Thus the 
primary balance in (4.2) yields a requirement that there be no net torques 
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on fluid elements. 

Thus writing T = T(z) + T'(x,z) 

J(T' + z. ^ z) = 0 


The solution is 



which is less than 1 for r > p and typical values. In principle this T' 

could be compared with the refractive patterns in fig. 6, for example. 

2 

Since T' == 1/r it will not satisfy thermal equilibrium demanded 

2 

by V T' = 0. Thus a weak circulation, constrained by the static stability 
will be set up. 

We imagine that the flow will be unstable when "gravity" is inverted, 
so that vertical excursions of cold fluid will be reinforced. This is 
essentially an inviscid criterion (valid because and a are large) 

a 

9T 

equivalent to the zero voltage condition that < 0* Looking at the 
coefficient of T^' in eqn. (4.2) we see that instability will be possible 
if 

. 22 ^ < 0 


where 9 is the polar angle. When 9 = - ir/2 instability is most likely. 
It is possible when 



a-6„ aV 
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* 

Let Pjjj = ra = r and setting = 2p we find that for our little experiment 

AV = 6 kv 

This compares well with the observed transition value of about 7 kv. 

Thus we have demonstrated that electrohydrodynamic forces do behave 
qualitatively as predicted. But more importantly we have shown that flows 
at large field strengths can be visualized easily and data concerning 
fluctuation frequencies, wavelengths, and transition points may be obtained 
from sequential photographs resembling fig. 7. 
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6. Calculations and experiment on spherical electroconvection. 

We have performed several calculations and experiments investigating 
the thermally driven flow in a hemispherical cap. That is, the spherical 
annulus of Figure 1 is divided into a hemisphere by a rigid insulating 
boundary at the equator. Motion in the upper hemisphere is studied. 

It is driven by cooling the inner surface relative to the outer. Thus 
T. < Tq. The fluid motion is now in response to both terrestrial gravity, 
and to electrodynamic forces. When AV= 0 there still will be motion 
because with this distribution of heating 

7T X 0 

In all runs * 126, $ = 2.36 and T, * 0. We have focused on the meri- 
r a 

dional circulation during axi symmetric (no variation with longitude) 
steady flow, so that the velocity vector has only two components 

w = radial velocity = w(r,(J>) 
and V = poleward velocity = v(r,4)) 

where is the 'latitude', zero at the equator. 

Our main purpose in doing these experiments was: 

a) to see how electrodynamic effects modify the 'Hadley' circulation driven 
by the sideways heating when AV= 0; 

b) to demonstrate that flow motions can be visualized by particle tracking 
at modest fields. 

The theoretical solutions of the governing non-linear Navier- 
Stokes equations (2.17 - 2.21) were obtained numerically by P.A. Gilman, 
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using a modified computer model originally designed to study solar con- 
vection (see Gilman, 1972). The results for fixed P^, T_, B depend only 

I a 

on R, and y* All results of this section will be presented in terms of 

a 

„ , gc. AT 
''a kv 

based on terrestrial gravity , the normal Rayleigh number, and AV, the 
applied ac voltage. 

The experimental apparatus has a brass inner spherical shell into 
which is pumped cool fluid from a thermostated bath. Thermal contrasts AT 
range from 1 to 10°C. The outer shell is glass coated on the inside 
with a transparent conductive coating (Indium oxide, by Perkin Elmer). 

The whole shell is immersed in an outer bath maintained at temperature T^. 
The working fluid is silicone oil (Dow Corning) laced with a small amount 
of merlescent pigment, a dielectric compound of small platelets that 
align themselves with the shear or with the electric potential surfaces. 
Time photographs are taken of particle streaks. From the photographs 
both flow structure and velocity data can be recovered. 

Figures 8-10 show the computed meridional circulation for AV = 0. 
It is basically a Hadley circulation of one cell with cold fluid sinking 
and hot fluid rising up the outer sphere in response to terrestrial 
gravity acting on hot and cold fluid. As R, increases velocities increase 

a 

and the flow becomes more of a boundary layer nature with a tendency 
towards stagnation near r = 1/2, d) = 20°. At R. = 11,000 the thermal 

a 

distribution is nearly conductive except near the equator where there 
is radial outflow (fig. 11). At large R^ a stable stratification is 

a 

developing in the interior (fig. 12) with strong radial thermal gradients 




Figure 8. Meridional circulation in a hemispherical shell with T. = 

a 

B = 2.36. These parameters are held fixed in all subsequent results 




Figure 9. Meridional circulation for AV 
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only near the walls. 

Figs. 13-14 show streak photographs taken by shining a narrow 
slit of light in towards the vertical axis and positioning the camera at 
90° to the slit. The computed and observed circulation patterns are 
similar. The vacant 'hole' of figure 10 is not really stagnant; the 
motion is just slower than the cutoff level required for an arrow. 

Quantitative comparisons between theory and experiment are shown 
in figures 15 and 16. The agreement is good. It is very hard to get 
good data in the high shear regions near the walls because small errors 
in obtaining the radial position of the particles lead to large errors 
in the predicted velocities. 

Since the agreement was so good for AV = 0 we then went on to 
study electrodynamic effects. However, a problem immediately arose. 

Fig. 17a and 17b show predicted meridional velocities for two cases. 

The 15,000 volt case is predicted to be indistinguishable from the zero 
voltage case. The gravitational Hadley circulation is dominant at these 
voltages. This is in part a consequence of having weak meridional thermal 
gradients (i.e. fig. 11). Since electrodynamic forces act via 

VT X VE^(r) 

and since is much less than (recall that VT x V(J>£ has a part 
ga Ip cos(})) our original scale analysis of electrodynamic effects vs. 
gravitational buoyancy effects was an overestimate, in this particular 
geometric configuration. This problem did not arise in the rod experiment 

VT X V(|) ^ = 0 


where 
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Figure 16. Comparison of theoretical and experimental (dots) values of v at <|» » 30®. 
R = 45,000. 

A 
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nor would it arise in the zero gravity laboratory. 

Figs. 18 show that the corresponding experimental circulations 
are nearly indistinguishable for aV = 0 to 16 kv. Note how nicely the 
particles line up with the potential surfaces (e.g., compare figure 18b 
to figure 18a where the particles are w.inking at the camera). At fields 
of 8 kv or less particle tracking can be done. This is a major result. 
However, if one looks closely at fig. 18c, one sees evidence of 'dancing' 
particles near the inner sphere at the pole. The particles pick up 
space charge and move in the field. This motion sets a voltage upper 
limit to the use of this technique. 

We have shown that flow measurements are possible at the fields 
needed in zero-g. We would also like a conclusive demonstration of the 
electrohydrodynamic effects in this geometry (to complement our previous 
demonstration and those of other workers mentioned in section 3). We 
have used the numerical model to suggest what we should look for. Fig. 

19-20 show that at large field strength the circulation becomes multi- 
cellular. Near the pole cold fluid rises in response to the fact that 

2 4 

the "gravity" has been reversed near the inner sphere (recall that E - 1/r 
so electrodynamic forces are biggest at smaller r). There are also several 
'plumes' (more obvious in fig. 20). These are thought to be electro- 
convective plumes, again with cold fluid rising, coming out of the 
boundary layer where "gravity" is inverted and being swept towards the 
equator by the main Hadley circulation. The time dependence can be seen 
by comparing figs. 21a and 21a. The 'plumes' are actually toroidal rolls 
moving equatorward. 

It is suggested that future work be aimed at measuring the critical 
value of AV required for the first appearance of these electroconvecti ve 
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instabilities. For these axi symmetric modes the critical voltage is 
about 45 kv. This is beyond our capability at M.I.T. (We do have a 50 kv 
transformer at Colorado University where this work will be continued.) 
However, in the presence of the radial shear of the main Hadley cell the 
most unstable modes should not be axisyimetric toroids but poloidal 
banana rolls (e.g., w* = w'(r,6), 6 is longitude). Being more unstable 
they should occur at lower AV. During the continuation of this work 
at Colorado University we will compute the critical value of AV for 
poloidal rolls, and compare with experiment. Obviously the voltage will 
still be high; we anticipate something like 30 kv. Thus no particles 
can be used. We will measure the refraction of a laser beam reflected 
off the inner sphere as the apparatus turns slowly (T, « 1). The rapid 

a 

fluctuation of the reflected beam will indicate the presence of non- 
axi symmetric rolls. The critical value of AV and the observed wavelength 
will be compared with theory, thus obtaining an accurate measure of 
electroconvective effects in this geometry and checking out another 
possible visualization scheme. 
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7. Conclusions and recommendations for zero-gravity experimental design. 

We have seen how terrestrial gravity effects swamp out electro- 

convective effects, and limit the application of this geophysical flow 

simulation method to low-g laboratories. Since the convective time 

scales are large an orbitting laboratory is preferred. It is easily 

seen that the apparatus cannot be too big because E- fields fall off as 
2 

1/r . Thus the optimum size for the spherical annulus is a few centi- 
meters for and one or two cm for d. But as we saw in section 4, 
the Taylor number T^ can still be large enough so that experiments can 

a 

be conducted in the rotational ly dominated regions of parameter space. 

Once the size (a few cm), temperature (~ 30°C) and voltage (< 8 kv) 

are set up, one might ask how big a feed-through neck is needed (see fig. 1) 

Alternatively one might ask how far away from the neck could one conduct 

experiments. The neck feeds high-voltage and coolent to the inner 

sphere so it must be on the order of a few millimeters in radius. Near 

2 

the neck there will be fringing fields. If E is not a function only 
of r there will be 'spurious' electroconvective buoyancy forces which 
are not radial, and hence that have no geophysical analogs. Thus we 
can only operate geophysical fluid dynamics experiments in that region 
of the spherical annulus where E is almost radial, or where the electric 
potential <(> is radial. 

We have tried to estimate the neck influence region by solving 
Laplace's equation for the geometry of fig. 22. An inner sphere is at 
potential (j> = 0, and an outer sphere with a hole in it at angle $ from 
the north pole, is at potential (j> = V = 1. Let 6 be the angle measure 
from the pole downwards to the hole. Clearly we will be best off if 
0 = $ -»- 7t/2 or d 0. Figure 23 shows computed <p fields for several 









-61- 


Fig. 23. Potential distributions. Horizontal axis is x = cos0. Neck is 
at X = 1, the left hand margin. Vertical axis is r. In a) Rq/R^ = 1.25, 
cos$ = -1. In b) Rq/R^ = 1.25, cos$ = -.5. In c) Rq/R^- = 1.25, cos$ = 0. 
In d) Rq/R^ = 2, cos$ = -.75. In e) Rq/R^- = 8, cost = -.75. xhe right 
hand axis Is at the north pole . Midway across each plot Is the equator. 
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cases. It is obvious that if we confine ourselves to experiments in a 
hemisphere that there will be no problem (our typical experiment has 
Ro/Ri - 1.3, cos$ == -.83). Looking at fig. 23b, equatorial annulus 
studies could conceivably be done to 45® S. 

We have explored several visualization schemes for obtaining 
data on simulated flows in the spherical annulus geometry with electric 
fields applied. Some typical results are apparent from earlier photo- 
graphs. They are: 

a) particle tracks or streaks: useful for meridional flow structure 
and velocity. Can be used for low velocities (.1 mm/sec or less as 
observed here in fig. 14). Disadvantage: need to be injected if g 0; 
otherwise will settle onto and adhere to surfaces. Only good for aV = 8 kv 
or less. 

b) Moir^ (refractive index) photos: useful for planforms, wave- 
lengths, frequencies and critical vlaues. Difficult for quantitative data. 
Good to breakdown point of the fluid. Best method for obtaining global 
structure of flow in longitude and latitude. Disadvantage: none. (Note 
other variants of this method may be better, e.g., shadow graph/Schlieren. 
Needs further technical study.) 

c) Remote anemometry (Laser doppler): useful for direct measurement 
of flow velocity at single point. We have used this method with the 
spherical annulus and found the method to be reasonably easy to set up and 
align and gives accurate results if velocities ~ 1 cm/sec are being measured. 
Otherwise vibration and alignment problems lead to noise problems. Consider- 
ing the care required and the complexity of such systems, it is not recom- 
mended for space use, although it can be a useful tool for exploring the 
fluid dynamics of flows in the terrestrial laboratory relevant to this 
problem. 



Figures 24 and 25 summarize current thought on a suitable zero-g 
experiment. It should be small and self-contained with the thermal 
'bath' for the inner sphere and the high voltage power supply located 
in the rotating housing. The cell with some rough dimensions is shown 
in fig. 25. We very much prefer the visualization method of b) above. 

It is simple, reliable, inexpensive, and will serve to give in a direct 
photographic record all the necessary data on fluctuation plan-form (wave 
axis tilt, for example), wavelengths, wave slopes, etc. The outer surface 
will be cooled or heated with an air jet. Different cells are suggested 
for mountain barrier problems, problems with non-uniform heating on 
spherical surfaces; our sketch shows the simplest configuration which 
undoubtedly should be the first to fly in order to test oiit under 
actual zero-g conditions the concepts we have outlined and studied in 
this report. 



Figure 24. Sketch of electroconvection apparatus. 



65 


glass hemisphere, transparent glass hemisphere, transparent 

conducting coating on outside conductive coating on Inside 



Figure 25. The electroconvectlcn cell. 
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